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Gravitational Waves from Numerical Mergers 

NR black-hole merger simulations produce waveforms decomposed into (spin-weighted) 
spherical harmonics: r^p 4 {t, r, 0, 0) = r) _2 Y^{0, 0). 

® We work with strain-rate h = f modal power oc /?|^ 

9 Each mode has an amplitude and complex phase: rhi^ = 

® A handful of modes dominate energy flux; mostly {i, ±£). 

® (2, ±2) is sufficient for detection] other modes are important for parameter estimation. 



Modal power for xi_uu (left), and X4_0 0 (right). 
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Dominant Frequency Behavior 

Baker et al. (2008) noted that many important modes have common rofaf/ona/ frequency 

^im = ^imt 



Rotational frequency for several {i, m) modes for 4:1 nonspinning merger. 
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Implicit Rotating Source Picture 


[Baker ef a/. (2008), Kelly ef a/. (2011)]: 

9 Most important WF modes had consistent 
rotational phases through 

merger. 

9 Best matches are for ^ = m modes. 

® Rotational frequency model is a smoothed 
“step function” to fundamental QNM 
frequency: 


fi(f) = 0 = fif(i - m, 


m = c 


^ ^ ae 




® Poor for times earlier than ~ 60M before 
merger. 



Fit of Q(f) for (2, 2) mode of xi_uu (top) 
and xi_DD (bottom). 
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Modeling IRS Parameters 

9 Assemble broad set of aligned-spin BHB 
merger configurations. 

® Collect free parameters {C, n, a} over all 
configurations. 

® Easier to model cuts along BH parameter 
directions . . . 

® Symmetric mass ratio 

rj = M 2 /{M^ + M 2 f < 0.25 

« “Total” spin J s (q^j^ +j 2 )l(q^ + 1 ) 

9 Simplest fit model is product of mass-ratio 
and spin forms: 

C(r/J) = g{v)-hO) 

giv) = go + g^ ivo -v) + gzim - vf, 
h{j) = -\+h^j + h2f. 




BHB configurations used, i 
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Modeling IRS Parameters 


9 Assemble broad set of aligned-spin BHB 
merger configurations. 

® Collect free parameters {C, n, a} over all 
configurations. 

® Easier to model cuts along BH parameter 
directions . . . 

® Symmetric mass ratio 

rj = M2/{M^ + M2f < 0.25 
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“Total” spin j = {q^j^ +k)/{q^ + 1 ) 

Simplest fit modei is product of mass-ratio 
and spin forms: 

C(??J) = 9iv)-h0) 

9iv) = 9o + 9^ ivo - v) + 92 (vo - vf , 



h(j) = ^+h^J + h2f. 


Fit to C [top] and k [bottom] over all 
configurations. 
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Modeling IRS Parameters 

Remaining IRS parameters depend on the end-state Kerr hole: af) & b{Mf, af). 

9 We want full initial-configuration prescription: 



{M ^ , M 2 , Si , S 2 } ^ {^2f, b} 
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Modeling IRS Parameters 

Remaining IRS parameters depend on the end-state Kerr hole: af) & b{Mf, af). 

9 We want full initial-configuration prescription: 



fAyf Ayf o O'! encl-statG model r., ^ QNM theory . . 

, M 2 , Si , S 2 } ^ {Mf, af} ^ {Qf, b} 

® Many prescriptions available, covering different ranges (nonspinning, aligned-spin, 

generic-spin), e.g. Lousto etal. (2010), Tichy & Marronetti (2008), Rezzolla etal. (2008), 
Barausse & Rezzolla (2009), Lousto & Zlochower (2014) . . . 

9 Use the simplest prescription consistent with aligned-spin BHBs: 

Mf =^ - rjEisco - E2rf- - E^rj^ 

2 

- (1 + < 7 ^ 71 ) + Es{^ - q){j2 ~qh) + Ea(J2 + Q/i f) + Ed(J2 ~ < 7 /i f 

^ =h = 3 + 3 v{sa 3 + ssv + h) + ri(2y^3) + tzn + 

® ... not necessarily the best prescription. 
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Modeling IRS Multimode Amplitudes 

IRS model [Kelly et al. (201 1 )] suggested a general form for IRS amplitude functions: 




Three free parameters for each (£, m) pair: Aq, , az- 



m 


1 + a-i — 

f4) + 32 - fi) 


^22 ^33 ^44 



^ 33 , and A 44 for 4:1 nonspinning merger. 
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Modeling IRS Amplitude Parameters 



Amplitude parameters don’t generally work well with quadratic forms; use leading-order 
post-Newtonian 77-scaling as guidance: 


rll22 OC TTf 


^ , 107 5577 . ^ 

1 ( rz — I 7^ I + O 3/2 


42 42 

r/733 oc 77 \/l — 477 + x ^/^(-4 + 277) + O2 


r/744 oc 77 


x (1 - 377) + 




22 


593 1 27377 , 2 » ^ 

- 1 7577^ + O3 


GOOD 

GOOD 

BAD 


^0,22 ^0,33 ^0,44 



7^0, 22= 7\o, 33, and ^0,44 parameters for nonspinning configurations. 
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Other problems with A44 

Something funny happens to A 44 near q = 2 ... 


A22 vs A 
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. . . fall-off too quick for n = 0 QNM (1 /r 0.08); too slow for n = 1 QNM (1 /r 0.24) — 
mode-mixing? 
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Dealing with the (3, 2) Mode 

Kelly & Baker (2013) showed that “observed” (3, 2) mode at merger is largely (2, 2) mode, leaked 
through mismatch between sp/7er/ca/ harmonics -2^^ and sp/7e/'o/c/a/ harmonics -23^^- 


® QNM eigenfunctions are spheroidal 
harmonics 

9 Overlap with spherical harmonics is 

Se'em = fan -2y?{afa22-. 0 , 4 >) -2 Y^{e, 4 >)* 
® . . . leading to mixing coefficients 

_ a ^/22 

Pbasis,£2 — 

S2/22 

® Fits observed (3, 2) modes very well; 

(4, 2) numerics too uncertain 

® Can/should project a// WFs onto spheroidal 
basis? 



Amplitude and frequency for several runs, 
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Dealing with the (3, 2) Mode 



Kelly & Baker (2013) showed that “observed” (3, 2) mode at merger is largely (2, 2) mode, leaked 
through mismatch between sp/7er/ca/ harmonics -2^^ and sp/7e/'o/c/a/ harmonics -23^^- 


® QNM eigenfunctions are spheroidal 
harmonics 

9 Overlap with spherical harmonics is 

Se'em = fan -2y?{afa22-. 0 , 4 >) -2 Y^{e, 4 >)* 
® . . . leading to mixing coefficients 

_ a ^/22 

Pbasis,£2 — 

S2/22 

® Fits observed (3, 2) modes very well; 

(4, 2) numerics too uncertain 

® Can/should project a// WFs onto spheroidal 
basis? 
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Summary 



® Need: higher-quality, higher-mode waveforms — NRAR? SpEC? 

® Need: treatment of important i ^ m modes (e.g. (2,1), (3,2)) 

9 Need: better behavior of pre-merger IRS segment — John Baker’s talk 
9 Need: better treatment of amplitudes in general 
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